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^ ; 1 Introduction 

-4— » 

The problem on possible dimensions of the maximal Lie invariance algebras of differential equa- 
tions from a fixed class has a long history. Already S. Lie obtained exhaustive results concerning 
maximal dimensions of such algebras for ordinary differential equations (ODEs) of any fixed 
^ . order |14l S. 294-301]. Namely, he first proved that any first-order ODE possesses an infinite- 

C^^ \ dimensional Lie invariance algebra, the dimension of the maximal Lie invariance algebra of any 

*yN ■ second-order ODE (resp. any mth order ODE for m ^ 3) is at most eight (resp. not greater than 

^^ . m + 4), and these bounds are exact. Later these results were repeatedly reinvented, see e.g. [S]. 

ff^ \ Analogous results for systems of ODEs are much less known. We discuss some of them, which 

O ■ are relevant to the subject of the present paper. Thus, according to the remarkable lecture notes 

by Markus [151 PP- 68-69, Theorem 44], any system of second-order ODEs 



x = f{t,x,x), (1) 

K> . 

Vh ■ where x{t) = (x^(i), . . . ,x'^(t))'^, x = dx/dt, x = dx/dt, possesses the maximal Lie invariance 

. . .' algebra of dimensions not greater than (n-|-2)^ — 1. This result was later reproved in [21 Sections 4 

and 5]. It was also shown therein that the maximal dimension (n -|- 2)^ — 1 = n^ -|- 4n + 3 is 
reached for systems reduced by point transformations to the simplest system 

x = 0. (2) 

The maximal Lie invariance algebra q^ of the system ([2]) is generated by the vector fields 

dt, dxa, tdt, x"'dt, tdx<^, x'^d^b, tx"'dt + x^x'^dx'^, t^ dt + tx^ dx<^ , 

and is isomorphic to the Lie algebra sl(n -|- 2, C) (resp. sl(n + 2,IR)) for the complex (resp. 
real) case, see e.g. [TT]. Here and in what follows the indices a, h, c run from 1 to n, i.e. 
a,b,c = 1, . . . , n, and we use the summation convention for repeated indices. It can be checked 



that the system ([2]) is invariant with respect to the general projective group of C""*" (resp. R"^^^ 
consisting of the point transformations [13' S. 554] 



X = T^ , z = U,...,n, 



an+ifix^ H h a„+i,„x" + a„+i,„+i 



where aooi oqi, ••., an+i,n+i are homogeneous group parameters and x^ = t, and g^ is the 
Lie algebra associated with this group. In fact, the number of essential group parameters is 
(n + 2)^ — 1 as, supposing a homogeneous group parameter nonzero, we can set it to be equal 
to 1, simultaneously dividing of the numerator and the denominator in the expression for each x* 
by this parameter and then assigning parameter ratios as new parameters. 

Fels proved O [6] that up to point equivalence the system ([2]) is a unique system of the 
form ([T|) which admits an (n^ + 4n + 3)-dimensional Lie invariance algebra, which was earlier 
known only for linear systems [11] . Recently linearization criteria for systems from the class ([1]) 
have been independently investigated in [HHKlT]. See also the references therein. The maximal 
dimension of the maximal Lie invariance algebras for normal systems of mth order ODEs was 
estimated in [HI [9l [101 for an arbitrary m ^ 3, and it is known for the case m = 3 that up to 
point equivalence the system x = is a unique system for which the dimension of the maximal 
Lie invariance algebra reaches the maximal value n^ + 3n + 3 for such systems [5] . 

In a series of recent papers the study of Lie symmetries of systems of n (n ^ 2) linear 
second-order ODEs with constant coefficients was recovered. Namely, the cases n = 2 and 
n = 3 were considered in ^22j. In [3j Campoamor-Stursberg corrected results of [22j (see also 
comments concerning [22] in [l6]), and studied the case n = 4 as well as systems associated with 
diagonal matrices without restrictions on n. Certain results on the dimensions of the maximal 
Lie invariance algebras of such systems in the case of arbitrary n and matrices of the general 
Jordan form were obtained in [4j. 

The study of symmetry properties of systems from the class ([1]) is required by numerous appli- 
cations in mechanics, gravity, etc. Unfortunately, there are no general results on Lie symmetries 
of these systems. This is why even linear systems with constant coefficients are good objects 
for a preliminary investigation in spite of the well-known simple algorithm for the construction 
of their general solutions. Group classification of linear systems with constant coefficients gives 
examples what Lie algebras of vector fields are admitted by systems from the class ([T]) as their 
maximal Lie invariance algebras and what dimensions of these algebras are possible. Note that 
the above knowledge is important for the problem on linearization of systems from the class ([1]) . 

The purpose of the present paper, which is inspired by the papers [3 JH [161 [22] , is to exhaus- 
tively describe Lie symmetries of systems of second-order linear ordinary differential equations 
with constant coefficients over both the complex and real fields. We essentially enhance and 
generalize the results of [31 [H [HI [22] using a simple but effective algebraic approach. In par- 
ticular, we explicitly describe the maximal Lie invariance algebras possessed by systems under 
consideration with no restrictions on the number of equations and the form of system coefficients 
and derive the exact lower and upper bounds for the dimensions of these algebras. 

2 Main result 

Following [31 m \W[ |22] , we consider systems of linear second-order ODEs of the normal form 

x = Ax + Bx + C{t) (3) 

over the complex field. Here C{t) is a smooth n-component vector-function of t, and A and B 
are constant complex- valued matrices of dimension n x n, n ^ 2. Note that the choice of the 
underlying field (C or M) is not principal. We choose the complex field in order to make the 
presentation clearer. 



It is commonly known (see, e.g., [TT]) that the change of dependent variables x = exp(2^i)l/+ 
a;p(t), where Xp{t) is a particular solution of ([3]), maps ([3]) into the system 

y = Dy, where D = B — A^. 

By J we denote the Jordan normal form of the matrix D. Then there exists a nondegenerate 
matrix P such that D = P^^JP, and the point transformation y = Pz reduces the last system 
to the form z = Jz. As a result, for the study of symmetry properties of normal systems of 
linear second-order ODEs with constant coefficients it suffices to consider only the systems of 
the form 



X 



Jx, 



where J is a Jordan matrix, 

s 
K^ "J Xl 1 ki -\- ■ ■ ■ -\- Ks 



J 



n, 



(4) 



(5) 



1=1 



J^ is the Jordan block of dimension ki with the eigenvalue A^, 



L A; J ij 




if j = i, 
if j-i = 1, 

otherwise. 



«,J 



:ki, 



i.e., it is ki x ki matrix of the form 
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(A — k)^^ 1 ■ ■ ■ , (A — Xs)'^" are elementary divisors of the matrix J. 

In what follows, we also use the notation diag(7i, ... ,7;) for the / x / diagonal matrix with 
the elements 71, . . . , 7/ on the diagonal, and E = diag(l, . . . , 1) denotes the I x I unit matrix. 
Subscripts of functions denote differentiation with respect to the corresponding variables. 

Remark 1. If the matrix J is proportional to the unit matrix, J E (^"), the associated 
system (jj]) is reduced to the form (l2|) by a point transformation, which is called generalized 
Arnold's transformation [11] . This is why we exclude such cases from the further consideration. 
At the same time, the representation (l5|) of J as a direct sum of Jordan blocks may contain 
partial direct sums of copies of the same 1x1 Jordan blocks, and such partial direct sums are 
proportional to the unit matrices of the corresponding sizes. 

To compute Lie symmetries of system ^ , we use the standard Lie approach [181 US] • Acting 
by the second prolongation of the vector field 



Q = C{t,x)dt + r]''{t,x)d^a 

on the system Q and substituting second-order derivatives using 
condition 



we derive the invariance 



-2{Ct + Cx^x1){Jx)' = {Jrj)'. 



(6t + 2CxHX^ + ix-x-x1x1 + i^a {Jxf)x\ 



(6) 



Splitting (^ with respect to the derivatives x" gives the system of determining equations for 
the coefficients ^{t,x) and r}'^{t,x), 

Cx'^x'^ = 0, (7) 

Vx^x'^ = 0, a^b^c, rila^b = ix-^t, a ^ b, rj^^t^b = I^^h^ (8) 

Tf',^^ = i^a{Jx)\ a^b, 2ii,^ = iu + 2iAJx? + ix'^{J^)\ (9) 

ri^ + vliJ^r - 'i^tiJxf = {Jv)\ (10) 

where there is no summation with respect to the repeated index b. Then equations dTj) and dH]) 
simultaneously imply 

e = r(i)x" + fit), rf = Q{t)x-x'' + ry^"(t)x'^ + rf^t), (11) 

where ^°', SP , r] and rj are smooth functions of t. Substituting the expressions (jlip into the 
equations ([9]) and (jlOp . and taking into the account the condition J ^ (E^), we obtain the more 
specific expressions 

where ci, cq, r/ are constants and rj are smooth functions of t. The detail derivation of these 
expressions is the following. The substitution of the expressions (jlip into equations Q and the 
subsequent partial splitting with respect to x result in the system 

where there is again no summation with respect to the repeated index b. Let S = {oj, i = 
1, . . . , m}, m ^ n, he the subset of integers from {1, . . . , n} for which ^"* ^ 0. For any fixed 
Oj € 5 we have the condition ^li H'^^ = {Jx) /x =: ^i = const whenever b ^ ai. If n > 2 and 
m ^ 2, the constants Hi, i = 1, . . . ,m, coincide to each other. If n = m = 2, the above condition 
is equivalent to the set of equations (Jx)^ = fiix^, (Jx)"^ = fi2x'^, Ctt — 1^2^^ and ^f^ = fJ-iC"^. 
Then equations ([12]) are reduced to the simple equations (/xi — H2)C = and hence again 
/ii = /X2 as ^^^'^ 7^ 0. This is why both the cases considered contradict the condition J ^ (-E"). 
If n ^ 2 and m = 1, only a single equation among (J12p is not satisfied identically, which, after an 
additional arrangement, takes the form (Ja;)"^ = fiix"'^ and implies, jointly with the equations 
{Jx)^ = fiix^ for b 7^ Oj, the condition J G {E"") contradicting the supposition J ^ {E"'). 
Therefore, the only possibility is m = 0, i.e., ^" = for any a. Differentiating determining 
equations (jlOp with respect to t and then splitting with respect to x, we in particular obtain 
the matrix equation ^^ttuE"' — 2^% J = implying ^^^ = in view of the linear independence 
of J and E"^. Then rfl"' = for any pair of the indices a and b. 
Further, from the determining equations (|10p we also have 

r,o=V, r?° = (r?^o,...,r?"r, 

i.e. 77'^ is an arbitrary solution of the system ([3|), and additionally we have the following matrix 
equation for the n x n matrix H = {rj ): 

HJ - 2itJ = JH. (13) 

If ^i = for any Lie symmetry operator of system 1^, equation P^ is reduced to the 
condition of commutation of H with J, JH = HJ. Therefore, we arrive at the Frobenius 
problem: define all matrices H that commute with a fixed matrix J. This is a standard problem 
in matrix theory, see e.g. [3 Chapter VIII]. 



If there exists a Lie symmetry operator Q of system ^ with ^t 7^ 0, the corresponding 
matrix H satisfies the inhomogeneous matrix equation (I13p 



JH -HJ = kJ, (14) 

where k = —2c\ ^ 0. This equation is compatible if and only if the matrix J is nilpotent, 
which directly follows from Lemma 4 of [12. p. 44] or from Theorem II of [21]. The matrix 
equation ()14p is an inhomogeneous linear system of algebraic equations with respect to the co- 
efficients of the matrix H. Hence the general solution of (|14p is represented as the sum of a 
particular solution of ()14p and the general solution of the associated homogeneous matrix equa- 
tion JH = HJ discussed in the previous paragraph. A particular solution of the equation (J14p is 
Kdiag(l, 2, . . . , fci, 1, 2, . . . , A;2, • • • , 1, 2, . . . , /cg), where /ci, k2, ■ ■ ■ , kg are the sizes of the Jordan 
blocks of the matrix J, cf. ^. 

Summing up, we obtain the following theorem. 

Theorem 1. Suppose that J is a matrix of the Jordan form ^ which is not proportional to 
the unit matrix. The maximal Lie invariance algebra q"^ of the system, x = Jx is 

(;f™, m = l,...,2n, 'H^,£=1,...,N, T), 
or 

(;f™, m= l,...,2n, n^,£ = l,...,N, T, V), 
if J is a non-nilpotent or nilpotent matrix, respectively. Here 

;f"* = ^™»(t)9^a, n^ = {H^f''x''d^,, T = dt, P = tO* - 27"V9^a, 

the vector- functions cp"^ = {ip"^^ (t) , . . . , ip"^^ (t)) , m = l,...,2n, form a fundamental set of 
solutions for the system x = Jx, H^ , i = 1, . . . ,N, exhaust linearly independent matrices that 
commute with the matrix J, and 7 = (7"^) = diag(l, 2, . . . ,ki,l,2, . . . ,k2, . . . , 1, 2, ... , kg). 

By N = N{D) we denote the number of linearly independent matrices that commute with 
an n X n matrix D. It is obvious that N{D) = N{D) if the matrices D and D are similar. 

Corollary 1. The dimension of the maximal Lie invariance algebra q^ of the system x = Dx 
with D ^ (-E") is 2n+A^+l or2n+N+2 if D is a non-nilpotent or nilpotent matrix, respectively. 

Let the matrix J of the form ([S]) be a Jordan form of the matrix D and aij denote the degree 
of the greatest common divisor of the polynomials (A — Aj) ' and (A — Aj) ^, i.e., aij = if 
Aj 7^ Xj and aij = min(A;j, kj) if Aj = Xj. Then the number N(D) is equal to [7, p. 221] 

s 
N=Y,^^J■ (15) 

Let Xi(A), ..., Tq{X) form the complete set of nonconstant invariant polynomials of the 
matrix D with the degrees ni ^ • • • ^ n^ > 0. Each invariant polynomial Xq,(A) is a product 
of co-prime elementary divisors, Xq,(A) = (A — Ai)'^"^ • • • (A — Xp)'^°'p, a = 1, . . . ,q. Here Ai, . . . , 
Xp are all the distinct eigenvalues of the matrix D, dij ^ d2j ^ • • • ^ dqj ^ 0, j = 1,2, ... ,p, 
Ua = dai + • • • + dap, OL = I,. . . ,q, ui + ■ ■ ■ + Ug = n, and, therefore, Xi(A) • • ■Tq{X) is the 
characteristic polynomial of the matrix D. According to [7, p. 222, Theorem 2], we have one 
more representation for A^ = N{D), 

N = ni+2,n2 + --- + {2q-l)nq. (16) 



We list other elementary properties of N(D), which are needed for the further consideration. 
Thus, N{D) = n mod 2, i.e. N{D) takes only odd (resp. even) values for odd (resp. even) n. 
The value of N{D) is completely defined by the tuple n = (ni, . . . ,nq) of the degrees of the 
nonconstant invariant polynomials of D or, equivalently, the partition of n into the integer 
summands ni ^ • • • ^ n^ > 0, n = ni + • • • + Ug. The representation (fT6]) implies that 
n ^ -/V(L') ^ n^ for any n x n matrix D. The equality N(D) = n holds if and only if gr = 1 and 
hence n = (n), i.e., all the elementary divisors of D are co-prime in pairs or, in other words, 
all the eigenvalues of D are distinct in pairs. The maximal value A^ = n^ is reached only if the 
matrix D is proportional to the unit matrix, D € (i?"), as then the number of (nonconstant) 
invariant polynomials associated with D is also maximal and equals n, and n = (!,...,!). The 
submaximal value of N(D) equals N = n'^ — 2n + 2, and it is attained only if n = (2, 1, ... , 1); 
then D is similar to either j|^ (©"Ji^ J^J or ( 0"^]^ JlJ J^^, where Ai / A2. The next 
less value A^ = n^ — 4n + 8 corresponds to the tuple n = (2, 2, 1, . . . , 1), n ^ 4. 

Theorem [1] and Corollary [1] jointly with (jlSp and ()16p give a much more effective algorithm for 
computing the dimensions of the maximal Lie invariance algebras of systems from the class ([3]) 
than those existing in the literature, cf. [H Proposition 4]. We also can explicitly construct 
a basis of such algebra for a fixed Jordan matrix. In particular, we directly derive series of 
simple estimates for dimensions of these algebras (cf . also O H] ) • 

Corollary 2. The maximal Lie invariance algebra of the system, x = Dx is of minimal dimen- 
sion 3n + 1 among systems of the form ^ if and only if the matrix D is not nilpotent and all 
the elementary divisors of D are co-prime in pairs. 

In other words, the Jordan matrix associated with D consists of either a single Jordan block 
with a nonzero eigenvalue or a few Jordan blocks with eigenvalues different in pairs. If D is 
similar to the single Jordan block Jq with the zero eigenvalue, the dimension of the maximal 
Lie invariance algebra of the system x = Dx equals 3n + 2. 

Corollary 3. The dimensions of the maximal Lie invariance algebras of the systems of the form 
X = Dx, where D ^ {E"'), are not greater than n^ + 4 and this upper bound is reached if and 
only if the matrix D is similar to the nilpotent Jordan matrix Jg ( ©"=1 Jq) ■ 

Therefore, for any system from the class ([3]) which is not similar to the simplest system x = 
the dimension of its maximal Lie invariance algebra 3 satisfies the inequality 

3n + 1 ^ dim g ^ n^ + 4, 

and the bounds are exact. The value dimg = n^ + 3 is attained for any system x = Dx with 
the matrix D similar to either Jx (B [ 0"=i Jx ) ; where Ai 7^ 0, or ( 0"j^i^ "^a ) ® "^A ' where 
Ai 7^ A2, and only for elements from the equivalence subclasses of such systems with respect to 
point transformations. As above mentioned, we have dimg = 3n-|-2 if D is similar to Jq. At the 
same time, only in the case 2 ^ n ^ 4 for each integer value p from the interval [3n + 1, n^ + 4] 
there exists a system from the class ([3]) whose maximal Lie invariance algebra is of dimension p. 
For n ^ 5, there exists, in particular, no n x n matrix D with dimg G [n^ — 2n + 11, n^ + 2]. 
The number of such exceptional intervals grows with increasing n. 

If the matrix J is diagonal, i.e. all its elementary divisors are of degree 1, and additionally 
J ^ (E^) then we have A^ = N{J) = J2^=i '''h where Ai, . . . , Ap are all the distinct eigenvalues 
of the matrix J and rj is the multiplicity of Aj, i = 1, . . . ,p. This is why Corollary [T] directly 
implies Proposition 3 of |3]. 

Given a Jordan matrix J with a single eigenvalue A (i.e. J = J^^ © J^^ © • • • © J^", where 
ki ^ k2 ^ ■ ■ ■ ^ ks and fci + • • • + /cg = n), it follows from (jlSp and ()16p that 

s s— 1 s 

In view of Corollary [H this essentially simplifies Theorem 2 of |3] . 



Remark 2. Given a matrix D with real entries, the number of hnearly independent solutions of 
the systems x = Dx (resp. the number A'^ = N(D) of linearly independent matrices commuting 
with D) over the real field is the same as over the complex field. Therefore, all results of this 
section are directly extended to the real case. Thus, in Theorem[T]it suffices to consider real Jor- 
dan matrices and to take real counterparts for the vector-functions ip"^ = {(p'^^{t), . . . , ip'^^(t)) , 
m = 1, . . . , 2n, and the matrices H^, i = 1, . . . ,N. See the second example in the next section. 



3 Illustrative examples 

We present two simple examples that illustrate Theorem [TJ 

Example 1. For system Q with the Jordan matrix J = J^ (B J^ , i.e., the system 



X = Xix + X 



2 ■■2 

X 



Aix^, x^ 



X2X + X , X = X2X , 



(17) 



there are two different cases subject to the eigenvalues Ai and A2, namely Ai / A2 and Ai = X2- 
For Ai 7^ A2 the single nonconstant invariant polynomial of the matrix J is (A — Ai)^(A — A2)^. 
Hence A^ = 4, and any matrix commuting with J has the form [71 Chapter VIII] 



H 



Finding a fundamental set of solutions of the system ()17p and using Theorem [U we obtain the 
13-dimensional maximal Lie invariance algebras of the system (|17p 



f," 


ri'' 





\ 







V^ 














7?33 


^34 


Vo 








ri'V 



X d^i +x 8^2, X d^i, X 8^3 +x d^i, x 9^3, dt) 

for Ai 7^ A2 7^ or 

{e'^\td^i+2^/xld^2), e-'^\td^i-2^/x'id^2), e^^^d^i, e'^^^d^i, t^d^3+6td^4, 
t 8^3 + 28^4, td^3, d^3, X d^i +x 8^2, x d^i, x 8^3 +x 8^4, x 8^3, 8t) 

for Ai 7^ A2 = 0. 

If Ai = A2, the nonconstant invariant polynomials of the matrix J are (A — Ai)^ and (A — Ai)^, 
and hence A^ = 8. A matrix H commutes with J if and only if it has the form 



H 



As a result, we construct the 17-dimensional maximal Lie invariance algebra 

(e^*(t5,i+2v^9,2), e-^*(t5,i-2v^5,2), e^^^a^i, e-^*a,i, 

e^^*(t5^3+ 2^X744), e-'^'\t8^3-2^/x'i8^4), e^^*a^3, e-^^*9^3, x'^d^i + x^8^2, 

X 8^1, X 8^1 +x 8^2, X 8^1, X 8^3 +x 8^4, x 8^3, x 8^3 +x 8^4, x d^s, dt) 



i„- 


^12 


^13 


,»1 





v'' 





77I3 


^31 


^32 


^33 


^34 


Vo 


^31 





ri'V 



if Ai 7^ or the 18-dimensional maximal Lie invariance algebra 

{t^ d.j.1 + 6td.^2 , t'^d^i+2d^2, td^i, d^i, t'^d^s +6td^4, t'^d^.a +28^4, td^s, d^a, 
x^ d^i + x'^ 8^2 , x'^d^i, x^d^i + x^d^2, x^d^i, x^ 8^3 + x'^ 8^4 , x'^d^a, x^d^a + x^d^4, 
x^d^i, dt, tdt — 2j;"'^02.i — Ax'^d^2 — 2x^dx3 — ^a^^c?^:^) 

if Ai = 0. 

Example 2. Over the real field, consider the system 

Xi = fiXi + 1^X2, X2 = —l^Xi + 11X2, Xs = X4, X4 =0, X5 = 0, (18) 

which is of the form ^ with the real- valued Jordan matrix J = B? Jq © Jj , where R 



"flU 






A matrix H commutes with J if and only if it has the form 

\ 



/^ii 



H 



r] 



12 : 











-v'' 


v'' 
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:^33 


^34 


^35 
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^33 






V 



0:0 r]^^ : rf^ j 



Then Theorem [T] implies that the maximal Lie invariance algebra of the system (jlSp is generated 
by the 18 vector fields 

z^e"* cos /3t9^i +e"*((a^ - /3^ - /i) cos /3t - 2a/3 sin /3t) 0^2 , 
z^e"* sin /3i9^.i + e'^* ((a^ - f - ^) sin /3t + 2a/3 cos /3t)5^2 , 
i^e""* cos /3ta^i +e""*((a^ - /3^ - /x) cos /3t + 2a/3 sin /3t) 5^2 , 
z^e""* sin /3t5^i + e""* ((a^ - /S^ - /i) sin /?t - 2a/3 cos /?t) 9^2 , 

X 0^1 + X 9^2, X 5^1 —X dr^2, X 0^3 + X 93.4, X 53,3, X S^jS, x 5^5, x ^^.s, 5t, 



where a = (/x^ + z^^)^/"^ cos (^ arctan ^) , /3 = (^^ _|_ ^^2^1/4 gj^^ ^1 Cretan ^) . 



4 Conclusion 

In this paper we exhaustively study Lie symmetries of systems of second-order linear ordinary 
differential equations with constant coefficients over both the complex and real fields. The 
explicit description of the maximal Lie invariance algebra of any system from the class ([3]) is 
presented by Theorem [TJ Corollary [1] jointly with formulas (|15p and ()16p give a simple and 
algorithmic tool in order for computing the dimensions of such algebras. In particular, we 
show that these dimensions are completely defined by the degrees of the nonconstant invariant 
polynomials of the corresponding Jordan matrices in the reduced form (jl]) . We also give estimates 
for possible values of these dimensions. In order to fix ideas, the consideration is carried out over 
the complex field but it can be directly extended to the real case with preserving all assertions, cf. 
Remark [5J The results obtained in the present paper amend and generalize those of [51 Bl [TCt [22] . 
The advantages and simplicity of the approach proposed are illustrated by examples. 

As a next step, we plan to study Lie symmetries of more general systems than Q, in par- 
ticular, classes of linear systems with variable coefficients. This will need the investigation of 



the equivalence groupoids of the above classes, the description of their maximal normalized 
subclasses and combining the algebraic and compatibility methods of group classification |20] . 
A still open problem is also the linearization of systems ([T|) by point transformations. 
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